Quantum radiation reaction effects in multiphoton Compton scattering 
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Radiation reaction effects in the interaction of an electron and a strong laser field are investigated in the 
realm of quantum electrodynamics. We identify quantum radiation reaction with the multiple photon recoils 
experienced by the laser-driven electron due to consecutive incoherent photon emissions. After determining 
a quantum radiation dominated regime, we demonstrate how in this regime quantum signatures of radiation 
reaction strongly affect multiphoton Compton scattering spectra and that they could be measurable in principle 
with presently available laser technology. 
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PACS numbers: 12.20.Ds, 41.60.-m 

Radiation reaction (RR) represents one of the oldest and 
still open problems in electrodynamics fljj, |2[]. In classical 
electrodynamics a charged particle in the presence of an ex- 
ternal electromagnetic field emits electromagnetic radiation, 
continuously looses energy and momentum so that the emis- 
sion modifies the trajectory of the particle itself. To take into 
account this "reaction" of the emitted radiation on the elec- 
tron motion one should solve Lorentz equation together with 
Maxwell equations self-consistently. This amounts in solv- 
ing the so-called Lorentz-Abraham-Dirac (LAD) equation of 
motion [2]. This equation, although solving in principle the 
problem of RR classically, is plagued by inconsistencies, like 
the existence of the so-called runaway solutions [1, 2]. How- 
ever, these inconsistencies are fictitious in the realm of clas- 
sical electrodynamics as they arise at length and time scales 
well in the quantum regime fljj] . Moreover, in classical elec- 
trodynamics the LAD equation can be consistently approxi- 
mated by the so-called Landau-Lifshitz (LL) equation, which 
does not show the shortcomings of the LAD equation ll|,|3|,|4]]. 
The LL equation has not yet been tested experimentally and 
proposals have been put forward to achieve this goal [5, 6]. 
An alternative classical equation of motion including the RR 
phenomenologically has been proposed in Ref. |7]. 

Though in the realm of classical electrodynamics the force 
on the electron due to RR is much smaller than the Lorentz 
force in the instantaneous rest frame of the electron, it can be 
the dominant force in the laboratory frame [l|]. Accordingly, 
the so-called (classical) radiation dominated regime (RDR) 
for Thomson scattering can be identified 15l la]. In the case 
of an external laser field, the RDR is entered when the en- 
ergy radiated during one driving laser period is of the order of 
or exceeds the initial electron energy. For a plane wave laser 
field with carrier angular frequency ujq and electric field am- 
plitude Eq, and for an electron (mass m and charge — e < 0), 
initially counterpropagating with the laser field, with initial 
energy e an d momentum p > 0, the average energy emit- 
ted by the electron per unit time is of the order of am 2 Xo, 
where a — e 2 « 1/137, xo = (^o + Po)Eo/mE cr , with 
E cr = m 2 /e = 1.3 x 10 16 V/cm being the critical field of 
quantum electrodynamics (QED) (units with h = c = 1 are 
used throughout) [90. Accordingly, the classical RDR is char- 
acterized by the condition 15j,|8[]: Re = «Xo£o ^ 1> where 



£o = eEo/mu>Q is the relativistic field parameter, assumed 
here to be much larger than unity (present laser systems al- 
ready allow for values of £o ~ 100 1 10]). 

The above considerations are classical. Quantum effects 
become important if Xo ZL 1- I n f act > me motion of an ul- 
trarelativistic electron in a laser field is quasiclassical ll ill . 
Quantum effects amount to the recoil experienced by the elec- 
tron in the photon emission which is of the order of Xo£o 1211 • 
However, quantum RR effects are not exhausted by single- 
photon recoil. In fact, quantum single-photon emission spec- 
tra at xo "C 1 turn into the corresponding classical spectra but 
without RR effects included. 

In the present Letter we investigate the quantum RDR in 
the interaction between an electron and an intense laser field. 
We identify quantum RR with the consecutive photon recoils 
in multiple incoherent single-photon emissions by the laser- 
driven electron, i.e. in the successive single-photon emissions 
each occurring in different radiation formation lengths. The 
quantum RDR is characterized by the two parametric condi- 
tions: Xo ^ 1 an d Rq = a £,o ^ 1- The first condition implies 
that the recoil in each photon emission is in general signifi- 
cant and it is fully accounted for by working in the framework 
of strong-field QED [9]. The second condition, instead, in- 
volves the quantum parameter Rq and implies that the aver- 
age number of photons emitted incoherently in one laser pe- 
riod can be larger than unity. We employ a new, microscopic 
approach to describe multiple incoherent photoemissions, al- 
ternative to the usual kinetic one Illl - ll3ll . which allows us 
to obtain emission spectra without solving the kinetic partial 
integro-differential equations. In this approach, the radiation 
process consists of multiple channels, each corresponding to a 
different number of incoherently emitted photons. The change 
of the electron state is taken into account consistently at each 
emission event which happens in a statistically uncorrected 
way. This is in contrast to the approach in [14] where, the sta- 
tistical character of the photon emission is neglected. We cal- 
culate numerically quantum emission spectra with and with- 
out RR and show how significantly the photon spectrum due to 
multiple emissions may differ from that predicted in the case 
of a single emission. This can represent a possible observable 
to measure RR effects in the quantum regime at optical laser 
intensities of the order of 10 22 W/cm 2 [ 10]. 



An ultrarelativistic electron with initial four-momentum 
p M = (s , 0, -p , 0) (s = y/m 2 + pi and p > 0) is con- 
sidered to collide head-on with a plane wave characterized 
by the four-vector potential Aq((J>) — (0,zAof(<f>)), where 
(f) = uj n (t — y), A = —E /uJo and /(</>) is an arbitrary func- 
tion with \f'(</>) | max = 1 (/' (<f>) = df{<t>)/d<t>). The external 
plane wave field can be taken into account exactly by quan- 
tizing the Dirac field in the Furry picture and by employing 
so-called Volkov states as electron states [10]. On the other 
hand, the interaction between the Dirac field and the photon 
field is taken into account perturbatively up to first order, i.e. 
in each formation length only the emission of one photon is 
taken into account. If Xo ;$ 1 [Xo ^ 1] the emission prob- 
ability of j photons in the same radiation formation length 
is approximately a J_1 [(aXo ) J ' _1 ] times the probability of 
emission of one photon [10], therefore we restrict ourselves 
here to xo "C 1/a 3 ^ 2 ~ 10 3 [16]. In general, it is conve- 
nient to label the particle state with initial four-momentum 
p^ = (s,p) = (\Jm 2 + p 2 ,p) via the quantity p- = e — p y , 
which is a constant of motion in the presence of the plane 
wave with four- vector potential Aq(4>). Also, we will always 
consider situations in which e ^ m £o> therefore the trans- 
verse degrees of freedom of the particles in the problem can be 
neglected and only probabilities integrated over the transverse 
momenta will be considered. We also consider probabilities 
averaged over the initial electron spin and summed over the 
final electron spin and the emitted photon polarization. As a 
starting point, we employ the expression for the probability of 
one photon emission dPpJ /dud<j> by an electron with initial 
four-momentum p^ in the external plane wave per unit of the 
laser phase <f> and per unit u = k_/(p_ — fc_), withu G [0, oo) 
1: 
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where fc_ = w — k y , with U L = (w,k) = (|fe|,fe) be- 
ing the four-momentum of the emitted photon, where K v (x) 
is the modified Bessel function of i/th order and where 
X {4>) = P -\E(<t>)\/mE cr , with E(4>) = E f(<t>) being 
the instantaneous electric field of the plane wave. Note 
that x(<t>) = (P-/Po-)xo|/'0)| = (p-/po,-)Xo{<t>), with 
Xo(4>) = Xo I /'(</>) I- The above expression of dPpJ /dud(f> is 
valid in the limit of a constant crossed field when £o ^ 1- 
However, since in this limit the radiation formation length of 
the photon production process is of the order of Ao/£o *C Ao, 
with Ao = 2ir/u>o being the central wavelength of the plane 
wave |9|], then the same expression of the probability can be 
employed for an external laser field varying in a space scale of 
the order of Ao , but with the instantaneous value of the electric 
field strength. In this respect, the instantaneous value p_ ((f)) 
at the moment of emission should be employed in Eq. (fl}, but 
in a plane wave the quantity p_ ((f)) is constant: p- ((f)) = p_. 
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FIG. 1: A diagrammatic visualization of the calculation of incoherent 
multiphoton emission by a laser-driven electron. The double electron 
lines indicate that the electron states are Volkov states and the dots 
represent the emission of more than two incoherent photons. At each 
photon emission step the electron incoming quantum number p_ is 
changed and determined by the previous step. 



For sufficiently long and/or strong pulses the quantity 



Pp _ — f^° d(f> J^° dudPp a _/dud(f) becomes larger than 
unity, and it cannot be interpreted as an emission probability. 
In this case multiple incoherent photon emissions take place, 
i.e. the electron emits many photons but each in a different 
radiation formation length. In the classical limit of emission 
of low-energy photons with respect to the initial energy 
of the electron, this apparent contradiction can be easily 
resolved following the procedure developed for the so-called 
"infrared catastrophe" (see 11611 ). The conclusion is that in 
the classical limit the quantity Pp _ is actually the average 
number of photons emitted (see also 1 13(1 ) and, therefore, 
it is not contradictory that it can be larger than unity 11611 . 
In the quantum regime we proceed analogously by taking 
into account, however, that at each photon emission step 
the electron incoming quantum numbers are changed and 
determined by the previous step (see Fig. 1). The total 
probability Pp _ that an electron with initial momentum 
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Pq emits N photons incoherently is given by P t 

f[du] N J[d(f>} N dP$_ 
where p } ■, = 

1,...,N - 1 and 
J °° du N f n °° dujv-i • 
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f_ d(f)\ has been introduced. 
Since Volkov states are normalized to unity, the proba- 
bility of no-photon emission is also unity at this stage. 
However, since the total probability that either no-photons 
or any number of photons are emitted must be unity, the 
actual total probability V Po _ of emission of N photons 
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• . Analogously, the 
probability of no-photon emission is actually l/Af Po _ . In the 
classical limit, i.e. when xo "C 1, then the main contribution 
to the integral in Pp _ comes from the region Uj <C 1 for 
j = 1,...,N, then the multidimensional integral in Pp J_ 
becomes factorisable and we obtain the Bloch and Nordsieck 
result: P^l = (P$_) N /N\ and N Po ._ = exp(P^_) 
11611 . In contrast to that, in the quantum regime the individual 
probabilities dPLL 1 _/dujd(j>j in Pp _ are not independent 
from each other, because, as we have mentioned, the initial 
electron momentum at each emission is different. The effects 
of RR arise here from the fact that the emission of each 



photon modifies the electron state and consequently the next 
emissions. Note that the quasiclassical motion of the electron 
between two emissions is taken into account consistently 
here, in the sense that the only dynamical quantity of the 
electron appearing in the equations, i.e. Pj,-, does not change 
between two emissions. Now, we introduce the quantum 
photon spectrum dS^/dzu including RR, of a photon 
with fc_ = njpo.- emitted by an electron with initial po,_ 
(for notational simplicity the index po - is n °t included in 
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In this expression the zth term in the sum over i corresponds 
to the photon with fc_ being emitted at fa. In general, 

dPp_ /dmdcj) = (p_/po,_)(p_/po,- _ -co)^ 2 dPfj /dudcj) 
and to = (p_/pQ t _)u/(l + u) < 1, while the step func- 
tion i?(x) ensures that the corresponding term vanishes if the 
electron has not enough energy to emit the photon with fc_ 
(uo = 0). Furthermore, pj t - — Pj-i,_ — tvpo,- lPj,- = 
Pj-i,-/(]- + Uj)] if the jth photon emitted is the on with fc_ = 
wpo,- [kj,- = pj_i,_Wi/(l + Uj)] for j = 1, . . . ,N - 1 
(po,- = Po,-)- Since we assumed that £o 3> m £o then 
po.- = &o — Po,y ~ 2eo an d the photons are mostly emit- 
ted along the negative y-direction, i.e. k_ = u — k y « 2w, 
then the quantity dS^ /dm is essentially the differential av- 
erage energy emitted in units of the initial electron energy. 
From the above discussion it is clear that in order to be in 
the quantum RR regime it must be: Pp _ > 1 (significant 
incoherent emission of many photons during the whole laser 
pulse) and \a ^ 1 (non-negligible photon recoil). Starting 
from Eq. (Q]) with p_ = po,_ and by performing the variable 

change v! = 2u/3xo(4 l ) m Pj>o,--> one sees that if xo ^ 1> 
then Pp ;_ ~ a(m 2 /woPo,-)xoA0 = a£oA0, where A<p is 
the phase interval corresponding to the pulse duration. Cor- 
respondingly, the analogous condition of being in the quan- 
tum RDR, i.e. Rq = a£o > 1, is obtained from the above 
one by setting A<fi ~ 1 (one period corresponds actually to 
Acj) = 2n). Therefore, when RR is not important, i.e. when 
a^oAcf) <C 1, the emitted spectrum is given by the quan- 
tity dS^ RR /dw defined as dS^ RR /dw = wdP$_/dw = 

w J [dfa i dPp _ jdwdfa . 

Below we compare the quantum spectra d<Sg R jdw and 
dSg RR jdvo also with the corresponding classical quantities. 
The expressions dS'™ RR /dvu and dS RR /dw of the classical 



spectrum without and with RR can be obtained starting from 
Eq. (fl} and from dSQ RR /dm for an electron with initial mo- 
mentum Pq in the limit £;_ <C po,_ (see also 111 ill ). While 
without RR it is p-{fa) = po,- = p n ° RR (4>), the analyti- 
cal solution of the LL equation found in [8] indicates that if 
RR is included, it is p-{fa) = p -/h(4>) = p RR (0), with 
h(fa) = 1 + (2/3)R c J^ OD dfaf 2 ((t)'). Therefore, we have 
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If one calculates classically the average en- 
ergy (w) emitted in one laser period as (w) s=s 
eo(27r)~ 1 J d(j) j„ dvjdS^ RR /dzud(t>, one obtains that 
the condition (w) > £o is fulfilled when ifc <^ 1> which is the 
classical RDR condition. In the classical regime, accounting 
for RR amounts essentially to subtract at each instant of time 
t + dt, from the four-momentum of the electron at time t 
the average four-momentum that it emits as radiation in the 
short interval dt. In contrast to that, in the quantum regime 
the loss of energy and momentum by the electron through 
the emission of a photon is intrinsically a probabilistic event. 
Only if the emission occurs, then the state of the electron and 
the successive emissions are modified. If one calculates (ui) 
starting from dS l Q RR /dru, one sees that quantum corrections 
are negative [9], thus the classical "averaged" treatment of RR 
is expected to overestimate the effects of RR when quantum 
effects become important. 

Below, we shall provide a numerical example in order to 
show that in principle the quantum RR regime can be al- 
ready experimentally realized with presently available laser 
technology. We consider the simple case of a two-cycle si- 
nusoidal pulse with /((/>) = — cos(0) and <j> £ [0, Air\, In 
this case the integrations in <f> in Eq. (O can be approximately 
performed analytically and the numerical calculation of the 
remaining integrals via the Monte Carlo method becomes 
feasible (numerical time-consuming Bessel functions can be 
avoided by employing an alternative equivalent representation 
of dPp a '_ /dudip as in Eq. (4.2) in [ 13]). Quantum and clas- 
sical spectra with and without RR shown in Fig. 2 clearly 
display different behaviors. Since classically the emitted fre- 
quencies do not have the meaning of photon energies, the clas- 
sical spectra allow for the unphysical emission at w > 1, 
also in the case that RR is taken into account. This is forbid- 
den quantum mechanically and the existence of the physical 
boundary w = 1 originates a typical quantum piling-up of 
the radiation by an ultra-relativistic electron towards w = 1. 
Inclusion of RR in the quantum regime has mainly three ef- 
fects: 1) increase of the spectral yield at low energies, 2) shift 
to lower energies of the maximum of the spectral yield and 
3) decrease of the spectral yield at high energies. The physi- 
cal reason is that, due to RR, the electron looses its energy by 
emitting several relatively low-energy photons and the proba- 
bility of emitting one photon in the region very close to w w 1 
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FIG. 2: (Color online) Multiphoton Compton spectra calculated 
quantum mechanically with (solid, black line) and without (long 
dashed, red line) RR. For the sake of comparison, the corresponding 
classical spectra with (short dashed, blue line) and without (dotted, 
magenta line) RR are also shown. The error bars in dSg^/dzu stem 
from uncertainties in the multidimensional integrations in Eq. l[2}. 
The inset shows a zoom of the spectral region where different curves 
cross. The numerical parameters are: eo = 1 GeV, ljo — 1.55 eV 
and To = Eq/8tt = 5 x 10 22 W/cm 2 , corresponding to £ = 154, 
Xo = 1-8 and Rq = 1.1. 



is less than if RR is neglected. Fig. 2 also shows that the clas- 
sical treatment of RR (short dashed, blue curve) artificially 
enhances the above three effects of RR and, as we have men- 
tioned, this is mainly due to the classical overestimation of 
the average energy emitted by the electron. Finally, in Fig. 2 
emission of up to sixteen photons is taken into account imply- 
ing the numerical integration of up to 15* h -fold integrals, the 
remaining higher-order terms having been estimated to give a 
contribution of about 2 %. 

Here we have not considered the exact space-time shape 
of the laser field. It is however clear from the above physi- 
cal considerations, that the laser's pulse shape will not change 
qualitatively our predictions. Concerning the observability of 
the discussed effect, one could expect that the emitted pho- 
tons, by interacting again with the laser field might create an 
electron-positron pair and start a cascade process 117]. How- 
ever, from Fig. 2 we see that the photon spectrum dSW/dvo 
has a peak at voq w 0.04. The probability of pair creation by a 
created photon with a given vo is roughly suppressed by a fac- 
tor r\{vo) = exp(— 8/3x 7 (n7)), with x^w) — voxo assumed 
to be much smaller than unity [9]. Since in the numerical ex- 
ample xo = 1-8, it is x 7 (o7o) = 0.073 and r/(-cuo) ~ 10~ 16 . 
Thus, we can conclude that most of the created photons, at 
not too large w's, will escape the laser pulse and reach the 



detector. Finally, we have calculated from Fig. 2 that about 
1.74 ± 0.17 and 1.00 photons are expected to be emitted at 
10~ 5 < vo < 10 -3 (corresponding to photon energies be- 
tween 10 KeV and 1 MeV) per electron with and without 
RR, respectively, amounting to a relative difference of about 
(43 ±11) % (note that for this example substantial devia- 
tions from Eq. (Q]i according to 111811 and consequently from 
the estimated photon yields only apply at vo < 10~ 8 ). This 
and the fact that laser intensities and electron energies con- 
sidered above are within the reach of present technology, al- 
low in principle for experimental measurement of the quantum 
RR effects with an all-optical experiment exploiting electron 
wake-field acceleration 1 10, 19]. 
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